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ABSTRACT 

We use cosmological N-body/gasdynamical simulations that include star formation and feedback to 
examine the proposal that scaling laws between the total luminosity, rotation speed, and angular momen- 
tum of disk galaxies reflect analogous correlations between the structural parameters of their surrounding 
dark matter halos. The numerical experiments follow the formation of galaxy-sized halos in two Cold 
Dark Matter dominated universes: the standard = 1 CDM scenario and the currently popular ACDM 
model. We find that the slope and scatter of the I-band Tully-Fisher relation are well reproduced in the 
simulations, although not, as proposed in recent work, as a result of the cosmological equivalence between 
halo mass and circular velocity: large systematic variations in the fraction of baryons that collapse to 
form galaxies and in the ratio between halo and disk circular velocities are observed in our numerical 
experiments. The Tully-Fisher slope and scatter are recovered in this model as a direct result of the 
dynamical response of the halo to the assembly of the luminous component of the galaxy. We conclude 
that models that neglect the self-gravity of the disk and its influence on the detailed structure of the 
halo cannot be used to derive meaningful estimates of the scatter or slope of the Tully-Fisher relation. 
Our models fail, however, to match the zero-point of the Tully-Fisher relation, as well as that of the 
relation linking disk rotation speed and angular momentum. These failures can be traced, respectively, 
to the excessive central concentration of dark halos formed in the Cold Dark Matter cosmogonies we 
explore and to the formation of galaxy disks as the final outcome of a sequence of merger events. Dis- 
appointingly, our feedback formulation, calibrated to reproduce the empirical correlations linking star 
formation rate and gas surface density established by Kennicutt, has little influence on these conclusions. 
Agreement between model and observations appears to demand substantial revision to the Cold Dark 
Matter scenario or to the manner in which baryons are thought to assemble and evolve into galaxies in 
hierarchical universes. 

Subject headings: cosmology: theory - galaxies: formation, evolution - methods: numerical 
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1. INTRODUCTION 

The structural parameters of dark matter halos formed 
in hierarchically clustering universes are tightly related 
through simple scaling laws that reflect the cosmological 
context of their formation. These correlations result from 
the approximately scalefree process of assembly of colli- 
sionless dark matter into collapsed, virialized systems. An- 
alytical studies and cosmological N-body simulations have 
been particularly successful at unraveling the relations be- 
tween halo mass, size, and angular momentum, as well as 
the dependence of these correlations on the cosmological 
parameters. The picture that emerges is encouraging in its 
simplicity and in its potential applicability to the origin of 
scaling laws relating the structural properties of galaxy 
systems (see, e.g., Dalcanton, Spergel & Summers 1997, 
Mo, Mao & White 1998 and references therein). 

One example is the relation between halo mass and size- 
a direct result of the finite age of the universe. The nature 
of this correlation and its dependence on the cosmologi- 
cal parameters is straightforward to compute using sim- 
ple spherical collapse models of "top-hat" density pertur- 
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bations that have been found to be in good agreement 
with the results of numerical experiments (see, e.g., Cole 
& Lacey 1996, Eke, Cole & Frenk 1996, Eke, Navarro & 
Frenk 1998, and references therein). 

A second example concerns the angular momentum of 
dark halos, which is also linked to halo mass and size 
through simple scaling arguments. Expressed in nondi- 
mensional form, the angular momentum of dark matter 
halos is approximately independent of mass, environment, 
and cosmological parameters, a remarkable result likely 
due to the scalefree properties of the early tidal torques 
between neighboring systems responsible for the spin of 
individual halos (Peebles 1969, White 1984, Barnes & Ef- 
stathiou 1988, Steinmetz & Bartelmann 1995). 

Finally, similarities in the halo formation process are 
also apparent in the internal structure of dark halos. A 
number of numerical studies have consistently shown that 
the shape of the density profiles of dark halos is approx- 
imately "universal"; i.e., it can be well approximated by 
a simple two-parameter function whose formulation is ap- 
proximately independent of mass, redshift, and cosmology 
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(Navarro, Frenk & White 1996, 1997, hereafter NFW96 
and NFW97, respectively; Cole & Lacey 1996; Tormen, 
Bouchet & White 1996, Huss, Jain & Steinmetz 1999a,b). 

How do these scaling properties relate to analogous cor- 
relations between structural parameters of disk galaxy sys- 
tems? This paper is third in a series where this ques- 
tion is addressed through direct numerical simulation of 
galaxy formation in cold dark matter (CDM) dominated 
universes. In spirit, these studies are similar to those of 
Evrard, Summers & Davis 1994, Tissera, Lambas & Abadi 
1997, Elizondo et al. 1999, although are conclusions differ 
in detail from those reached by them. The first paper in 
our series (Steinmetz & Navarro 1999, hereafter SN1) ex- 
amined the origin of the Tully-Fishcr relation under the 
assumption that star formation is dictated by the rate at 
which gas cools and collapses within dark halos. We were 
able to show that the velocity scaling of luminosity and 
angular momentum in spiral galaxies arise naturally in hi- 
erarchical galaxy formation models. 

Large discrepancies, however, were observed in the zero 
point of these correlations, a result we ascribed to the 
early dissipative collapse of gas into the progenitor dark 
matter halos and to the subsequent assembly of the final 
system through a sequence of mergers (Navarro & Benz 
1991, Navarro, Frenk & White 1995, Navarro & Steinmetz 
1997, hereafter NS97). We concluded then, in agreement 
with a number of previous studies, that agreement be- 
tween models and observations requires a large injection 
of energy (presumably "feedback" energy from evolving 
stars and supernovae) in order to prevent much of the gas 
from cooling and condensing into (proto)galaxies at early 
times, shifting the bulk of star formation to later times 
and alleviating the angular momentum losses associated 
with major mergers (White & Rees 1978, White & Frenk 
1991, Kauffmann, White & Guiderdoni 1993, Cole et al. 
1994). 

This hypothesis remains to date the most attractive 
path to resolve the "disk angular momentum problem". 
On the other hand, further analysis has shown that the 
zero-point offset between the model and observed Tully- 
Fishcr relations is actually due to the high central concen- 
trations of dark matter halos formed in currently popular 
versions of the Cold Dark Matter cosmogony, and there- 
fore is unlikely to be reconciled through feedback effects 
(Navarro & Steinmetz 2000, hereafter NS2). We inves- 
tigate these issues in detail in the present paper, which 
represents our first attempt at implementing a realistic nu- 
merical formulation of feedback within a full-fledged sim- 
ulation of the formation of galaxies in CDM halos, and at 
gauging its effects on the origin of disk galaxy scaling laws. 

The paper is organized as follows. Section 2 motivates 
our approach by reviewing the scaling laws linking the 
properties of dark matter halos and by comparing them 
with observational correlations. A brief description of the 
numerical procedure is included in §3; full details of the 
star formation, feedback prescription, and relevant tests 
will be presented in Steinmetz & Navarro (2000, hereafter 
SN3). Section 4 discusses our main results regarding the 
slope, scatter and zero-point of the Tully-Fisher relation 
and of the relation between disk angular momentum and 
rotation speed. We summarize and discuss our main rc- 

3 We use the term "density contrast" to refer to densities expressed in 
H(z) is the value of Hubble's constant at redshift z. We parameterize the 



suits in §5, and conclude in §6. 

2. DARK HALO AND DISK GALAXY SCALING LAWS 
2.1. Mass, Radius, and Circular Velocity 

The radial distribution of mass in dark halos has no ob- 
vious edge, so defining halo "sizes" is a somewhat arbitrary 
procedure. A plausible and useful choice is to associate the 
radius of a halo with the distance from the center at which 
mass shells are infalling for the first time. This "virial 
radius" is easily estimated from N-body simulations and 
imposes, by construction, a firm upper limit to the mass 
of the galaxy embedded inside each halo: baryons beyond 
the virial radius have not had time yet to reach the center 
of the halo and therefore cannot have been accreted by the 
central galaxy. 

Numerical experiments show that virial radii depend 
sensitively on the enclosed mass of the system, in a way 
consistent with the predictions of the spherical top-hat col- 
lapse model. The "edge" of a halo occurs approximately 
where the mean inner density contrast, A, is of order a few 
hundred. 3 This result seems to apply equally well to ha- 
los of all masses, and depends only weakly on cosmology 
through the density parameter, f2, and the cosmological 
constant, A, (Eke et al. 1996, 1998), 
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Fig. 1. — The I-band Tully-Fisher relation compared with the 
results of the numerical simulations. Dots correspond to the ob- 
servational samples of Mathewson, Ford & Buchhorn, (1992), 
Giovanelli et al (1997), Willick et al (1997), and Han & Mould 
(1992). Error bars in the simulated magnitudes correspond to 
adopting a Salpeter or a Scalo IMF. 

One consequence of this definition is that all halos iden- 
tified at a given time have similar densities, from where it 
follows on dimensional grounds that the total mass of the 
halo and its circular velocity at the virial radius must be 
tightly related. In terms of the circular velocity, Va, at 
the virial radius, ta, halo masses are given by 
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This power-law dependence on velocity is similar to that 
of the I-band Tully-Fisher relation linking the luminosity 
and rotation speed of late- type spirals, 
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(see solid line in Figure 1) a coincidence that suggests a 
direct cosmological origin for this scaling law. Eqs. 2 and 
3 can be combined to read, 
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As indicated by eq. 5, reproducing the observed slope 
of the Tully-Fisher relation entails a delicate balance be- 
tween /mdsk, Tj, and the ratio V^/Va- The simplest 
possibility is that the three parameters are approximately 
constant in all halos. This is the case argued by Mo et al. 
(1998), who suggest that / m dsk ~ 5x 10~ 2 , Tj s» 1.7 h, and 
Vrot/VA ~ 1-5 are needed in order to reproduce observa- 
tions of galaxy disks. Although plausible, this assumption 
is at odds with the results of the numerical experiments we 
present below (§3), so it is worthwhile considering a second 
possibility: that / m dsk, Tj, and V^/Va are not constant 
from halo to halo but that their variations are strongly 
correlated, in the manner prescribed by eq. 5. We explore 
now how such correlation may emerge as a result of the 
dynamical response of the dark halo to the assembly of a 
disk galaxy at its center. 



where we have introduced the parameter Mdisk to repre- 
sent the total mass associated with the galaxy disk and 
Tj = Mdisk/ Li is the disk mass-to-light ratio in solar 
units 4 

Combining eqs. 2, 3, and 4, we find that the fraction of 
the total mass of the system in the galaxy disk is, at z — 0, 
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This result reemphasizes our implicit assumption that 
within the "virial radius" galaxy systems are dominated 
by dark matter. Further insight can be gained by compar- 
ing Mdisk to the total baryonic mass within ta, M^^k = 
(fi& /QqJMa- Assuming that the baryon density parameter 
is ilb ~ 0.0125 /i~ 2 , as suggested by Big Bang nucleosyn- 
thesis studies of the primordial abundance of the light ele- 
ments (Schramm & Turner 1997), the fraction of baryons 
transformed into stars in disk galaxies is given by, 
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Fig. 2. — The disk mass fraction versus the ratio between disk 
rotation speed and halo circular velocity. The thick dashed and 
solid lines correspond to the constraint imposed on these two 
quantities by the Tully-Fisher relation (eq. 5) in the ACDM 
and SCDM scenarios, respectively. Dotted lines correspond to 
the relation expected for galaxies assembled in NFW halos of 
constant "concentration" parameter, as labeled. Constant disk 
mass-to-light ratios are assumed throughout; Tz = 2 in the 
upper panel and Tj = 1 in the lower one, respectively. 

In the interest of simplicity, we shall restrict our anal- 
ysis to a case where the disk mass-to- light ratio, Tj, is 
assumed to be constant, and we shall use the "adiabatic 
contraction" approximation to compute the dependence of 
the velocity ratio K-ot/^A on / m dsk- The resulting relation 
is sensitive to the detailed mass profile of the dark halo 
and to the assumed structure of the disk. Assuming that 
the model proposed by Navarro, Frenk & White (NFW96, 
NFW97) 5 is a reasonable approximation to the structure 
of the halo, and adopting exponential disk models with 
radial scale lengths consistent with the assumption that 
halo and disk share the same specific angular momentum, 
it is possible to derive analytically the dependence of the 
velocity ratio on / m dsk (see Mo et al 1998 for details). We 
show the results in Figure 2. The thick lines (solid and 
dashed) in this figure correspond to the constraint enun- 
ciated in eq. 5 for two different cosmological models (the 
"standard" Cold Dark Matter model, sCDM: with param- 
eters 57 = 1 and h = 0.5, and a low-density, flat ACDM 

4 Note that this definition of the disk mass-to-light ratio includes stellar and gaseous mass. However, most of the baryons in the Tully-Fisher 
disks we consider here are actually in stars (this is true for observed as well as for simulated galaxies) so that the stellar and "baryonic" disk 
mass-to-light ratio differ very little. We shall not discriminate between them throughout this paper. Further, we assume Mj(0) = 4.15 for all 
numerical values quoted in this paper. 

5 According to these authors, the density profile of a dark halo is well approximated by p(r) = S c pcrit ( r / r 's) _1 (l + r / r s)~ 2 , where r s is 
a scale radius and <5 C is a characteristic density contrast. For halos of given mass, Ma, this formula has a single free parameter, which can 
be expressed as the "concentration", c = r&/r a . The characteristic density contrast and the concentration are related by the simple formula, 
<5 C = (A/3) c 3 /(ln (1 + c) - c/(l + c)). 



Because by definition baryons outside the virial radius 
have yet to reach the galaxy, M^JjJ is a firm upper bound 
to the baryonic fraction transformed into stars, implying 
that /bdsk<l (White et al. 1993). 

The slope and zero point of the Tully-Fisher relation 
therefore implies that the fraction of the total mass (and 
of baryons) transformed into stars is a sensitive function 
of the cosmological parameters (through the product fi h, 
and A), of the stellar mass-to-light ratio, and of the ratio 
between the rotation speed of the disk and the circular 
velocity of the surrounding halo. 
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Constraints from the slope of the Tully-Fisher 
relation 
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model, with ft = 0.3, A = 0.7, and h = 0.7). The up- 
per and lower panels adopt two different values for the 
disk mass-to-light ratio, Tj = 2, and 1, respectively. The 
rightmost point in each of the thick lines corresponds to 
the maximum disk mass fraction allowed by the baryonic 
content of the halo, i.e., f™^ k = Qfc/Qo- Dotted lines in 
this figure are the "adiabatic contraction approximation" 
predictions for different values of the NFW concentration 
parameter, c. 

Figure 2 illustrates a few interesting points. The first 
one is that the disk mass-to-light ratio and the cosmologi- 
cal parameters determine in practice the range of halo con- 
centrations that are consistent with the zero-point of the 
Tully Fisher relation. Halos formed in the sCDM scenario 
must have c < 3 (5) if Yj ss 2 (1). This effectively rules out 
the sCDM scenario, since N-body simulations show that 
halos formed in this cosmology have much higher concen- 
trations, typically c ~ 15-20 (NFW96). A similar con- 
clusion was reached by van den Bosch (1999), who finds 
through similar considerations that the large V I0 t/V~A val- 
ues expected in the sCDM scenario effectively rule out this 
cosmogony. 

The low-density ACDM model fares better, because the 
higher value of h and the lower value of A in this model 
imply smaller / m dsk at a given value of the velocity ratio 
K-ot/^A- However, for Yj = 2 (1), concentrations lower 
than about ~ 5 (12) are needed. As discussed by NS2, 
high-resolution N-body simulations of halo formation in 
the ACDM scenario yields concentrations of order ~ 20, in 
disagreement with these constraints unless Yj <C 1. Con- 
centrations as high as this are similar to those found for 
sCDM, and are systematically higher than the values pre- 
dicted by the approximate formula proposed by NFW97. 
This explains the apparent disagreement between the con- 
clusions of NS2 and those of semi-analytic models that 
claim reasonable agreement with the observed Tully-Fisher 
relation (Mo et al 1998, van den Bosch 1999): the discrep- 
ancy can be fully traced to the lower halo central concen- 
trations and lower disk mass-to-light ratios adopted by the 
latter authors. If concentrations are truly as high as re- 
ported by NS2, agreement with the Tully Fisher relation 
require Y/ <C 1, in disagreement with estimates based 
on broad-band colors of Tully-Fisher disks (which suggest 
Tj » 2) and with mass-to-light ratio estimates of the solar 
neighborhood (see NS2 for details). 

The second important point that emerges from Figure 
2 is that, in order to match the Tully-Fisher relation, halo 
concentrations must be an increasing function of / m dsk (as- 
suming Y/ ^constant). Indeed, as / m( jsk increases the 
thick solid and dashed lines cross (dotted) curves of in- 
creasing c. Since, according to NFW96 and NFW97, con- 
centration depends directly on halo mass — low mass halos 
are systematically more centrally concentrated as a result 
of earlier collapse times — this is equivalent to requiring 
/mdsk to be a function of halo mass. This is actually con- 
sistent with simple disk formation models where the mass 
of the disk is determined by gas cooling radiatively inside 
a dark halo, as first proposed by White & Rees (1978), and 
later worked out in detail by White & Frenk (1991). These 
authors show that, if disks form by gas cooling within an 
approximately isothermal halo, disk masses are expected 
to be roughly proportional to V^/ 2 (White 1996), imply- 



ing / m dsk « V^ 3 ^ 2 oc M^ 1 ^ 2 . We show below (§4.1) that, 
although this dependence is stronger than found in our nu- 
merical experiments, the overall trend predicted is nicely 
reproduced in our numerical experiments. 

Finally, Figure 2 illustrates that the structure and dy- 
namical response of the halo to the assembly of the disk 
may be responsible for the small scatter in the Tully- 
Fisher relation. For illustration, consider two halos of 
the same mass, and therefore approximately similar con- 
centration, where the fraction of baryons collected into 
the central galaxy, / m dsk, differs substantially. Provided 
that /mdsk >0.02, where the "adiabatic contraction" dot- 
ted curves are approximately parallel to the observational 
constraint delineated by the thick lines, these two galaxies 
will lie approximately along the same Tully-Fisher rela- 
tion. Even if the concentration of the two halos were to 
differ greatly its effect on the scatter of the Tully-Fisher 
relation would be relatively minor: at fixed / m dsk, Kot/^A 
changes by only about 20% when c changes by a factor of 
two. 

To summarize, assuming that the disk mass-to-light ra- 
tio is approximately constant for Tully-Fisher disks, the 
slope of the I-band Tully-Fisher relation may result from 
the combination of three effects: (i) the reduced efficiency 
of cooling in more massive halos, (ii) the mass dependence 
of halo concentrations, and (iii) the dynamical response of 
the halo to the disk assembly. This also offers a natural 
explanation for the small scatter in the observed Tully- 
Fisher relation. We shall use our numerical simulations to 
test the verisimilitude of this speculation below. 

2.2. Circular Velocity and Angular Momentum 

Another similarity between the properties of dark halos 
and galaxy disks concerns their angular momentum. N- 
body simulations show that, in terms of the dimensionless 
parameter, A = J|£ , | 1 / 2 /GA/^ 2 , the distribution of halo 
angular momenta is approximately independent of mass, 
redshift, and cosmological parameters, and may be ap- 
proximated by a log-normal distribution peaked at around 
A ~ 0.05 (Cole & Lacey 1996 and references therein). (J 
and E are the total angular momentum and binding en- 
ergy of the halo, respectively.) 
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Fig. 3. — Specific angular momentum vs circular velocity of 
model galaxies compared with observational data. Data cor- 
respond to the samples of Courteau (1997), Mathewson et al 
(1992), and the compilation of Navarro (1999). Specific angu- 
lar momenta are computed from disk scalelengths and rotation 
speeds, assuming an exponential disk model with a flat rotation 
curve. 

The binding energy depends on the internal structure of 
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the halos but the structural similarity between dark halos 
established by NFW96 and NFW97 implies that E is to 
good approximation roughly proportional to MaV^, with 
a very weak dependence on the characteristic density of 
the halo. The specific angular momentum of the halo then 
may be written as (see Mo et al 1998 for further details), 



A 



JA 



A 1 / 2 H(z) 



2.8 x 10 3 



200 km s" 



H 

H(z) 



km s 1 h x kpc, 



A 

200 



-1/2 



(7) 



where we have used the most probable value of A = 0.05 
in the second equality (see dotted line in Figure 3). The 
simple velocity-squared scaling of this relation is identical 
to that illustrated in Figure 3 between the specific angular 
momentum of disks and their rotation speed, 



j disk « 1.3 x 10 3 



200 km s _1 



km s _1 /i _1 kpc (8) 



(solid line in Figure 3), suggestive, as in the case of the 
Tully-Fisher relation, of a cosmological origin for this scal- 
ing law. 

Combining eqs. 7 and 8, we can express the ratio be- 
tween disk and halo specific angular momenta at z = 

as, 

1/2 /Tr \ 2 
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If the rotation speeds of galaxy disks are approximately 
the same as the circular velocity of their surrounding ha- 
los, then disks must have retained about one-half of the 
available angular momentum during their assembly. 

The velocity ratio may be eliminated using eq. 6 to ob- 
tain a relation between the fraction of baryons assembled 
into the disk and the angular momentum ratio, 
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This combined constraint posed by the Tully-Fisher and 
the angular momentum-velocity relation is shown in Fig- 
ure 4 for two different cosmological models. As in Figure 
2, thick solid lines correspond to the "standard" cold dark 
matter model, sCDM, and thick dashed lines to the ACDM 
model. Each curve is labelled by the value adopted for the 
disk mass-to- light ratio, T/. The precise values of /bdsk 
and fj along each curve are determined by and the ratio 
K-ot/^A, and are shown by starred symbols for the case 
V m t = V A and Tj = 1. 




Fig. 4. — The fraction of the baryons assembled into a disk 
galaxy (/bdsk) versus the ratio between the specific angular mo- 
menta of the disk and its surrounding halo {fj). Thick solid 
and dashed lines correspond to the constraints imposed by the 
Tully-Fisher relation (Figure 1) and by the relation between 
rotation speed and angular momentum (Figure 3). The solid 
(dashed) thick line corresponds to the sCDM (ACDM) scenario, 
shown for different values of Tj, as labeled. Symbols corre- 
spond to simulated galaxy models as per the labels in Figure 
1. 

One important point illustrated by Figure 4 is that disk 
galaxies formed in a low-density universe, such as ACDM, 
need only accrete a small fraction of the total baryonic 
mass to match the zero-point of the Tully-Fisher relation, 
but must draw a much larger fraction of the available an- 
gular momentum to be consistent with the spins of spiral 
galaxies. For example, if V Iot = Va and Tj = 1, disk 
masses amount to only about 30% of the total baryonic 
mass of the halo but contain about 60% of the available 
angular momentum. This is intriguing and, at face value, 
counterintuitive. Angular momentum is typically concen- 
trated in the outer regions of the system (see, e.g., Figure 
9 in NS97), presumably the ones least likely to cool and 
be accreted into the disk, so it is puzzling that galaxies 
manage to tap a large fraction of the available angular mo- 
mentum whilst collecting a small fraction of the total mass. 
The simulations in NS97, which include the presence of a 
strong photo-ionizing UV background, illustrate exactly 
this dilemma; the UV background suppresses the cooling 
of late-infalling, low-density, high-angular momentum gas 
and reduces the angular momentum of cold gaseous disks 
assembled at the center of dark matter halos. 

The situation is less severe in high-density universes 
such as sCDM; we see from Figure 4 that disks are re- 
quired to collect similar fractions of mass and of angular 
momentum in order to match simultaneously the Tully- 
Fisher and the spin-velocity relations. As a result, any 
difficulty matching the angular momentum of disk galaxies 
in sCDM will become only worse in a low-density ACDM 
universe. 

Note as well that the problem becomes more severe the 
lower the mass-to-light ratio of the disk. Indeed, from the 
point of view of this constraint, it would be desirable for 
disks formed in the ACDM scenario to have T/>2; in this 
case fj ~ /bdsk would be consistent with the constraint 
posed by observed scaling laws. However, this is the oppo- 
site of what was required to reconcile highly-concentrated 
halos with the zero-point of the Tully-Fisher relation. This 
conundrum illustrates the fact that accounting simultane- 
ously for the mass and angular momentum of disk galax- 
ies represents a serious challenge to hierarchical models of 
galaxy formation. 

3. NUMERICAL EXPERIMENTS 

3.1. The Code 

The simulations were performed using GRAPESPH, 
a code that combines the hardware N-body integrator 
GRAPE with the Smooth Particle Hydrodynamics tech- 
nique (Steinmetz 1996). GRAPESPH is fully Lagrangian 
and optimally suited to study the formation of highly non- 
linear systems in a cosmological context. The version used 
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here includes the self-gravity of gas, stars, and dark matter 
components, a three-dimensional treatment of the hydro- 
dynamics of the gas, Compton and radiative cooling, the 
effects of a photo-ionizing UV background (NS97), and a 
simple recipe for transforming gas into stars. 

3.2. The Star Formation and Feedback Algorithm 

The numerical recipe for star formation, feedback and 
metal enrichment is similar to that described in Steinmetz 
& Muller (1994, 1995, see also Katz 1992, Navarro & White 
1993). Full details on our present implementation are pre- 
sented in Steinmetz & Navarro (2000), together with vali- 
dating and calibrating tests. A brief description follows. 

"Star particles" are created in collapsing regions 
that arc locally Jeans unstable at a rate controlled 
by the local cooling and dynamical timescales, p* = 
^Pgas/meLxiTcooijTdyn). The proportionality parameter, 
c*, effectively controls the depletion timescale of gas, which 
in high-density regions, where most star formation takes 
place and where T dyn » t coo; , is of order Td yn /c* « 
(4 7 rGp 9as )- 1 /2 c -i. ' 

The equations of motion of star particles are only af- 
fected by gravitational forces, but newly formed stars de- 
volve ~ 10 ergs (per solar mass of stars formed) about 
10 7 yrs after their formation. This energy input, a crude 
approximation to the energetic feedback from evolving 
massive stars and supernovae, is largely invested in raising 
the temperature of the surrounding gas. As discussed by 
Katz (1992) and Navarro & White (1993), this form of en- 
ergetic feedback is rather inefficient; the high densities typ- 
ical of star forming regions imply short cooling timescales 
that minimize the hydrodynamical effects of the feedback 
energy input. The net result is that the star formation his- 
tory of an object simulated using this "minimal feedback" 
formulation traces closely the rate at which gas cools and 
collapses within dark matter halos (SN1). 

We have generalized this formulation by assuming that 
certain fraction of the available energy, e v , is invested 
in modifying the kinetic energy of the surrounding gas. 
These motions can still be dissipated through shocks but 
on longer timescales, leading to overall reduced star for- 
mation efficiencies and longer effective timescales for the 
conversion of gas into stars. 

We have determined plausible values for the two pa- 
rameters, c* and e v , by comparing the star forming prop- 
erties of isolated disk galaxy models with the empirical 
"Schmidt-law" correlations between star formation rate 
and gas surface density reported by Kennicutt (1998, see 
SN3 for full details). In the case of "minimal feedback" 
(e v = 0) low values of c*, typically ~ 0.05, arc needed to 
prevent the rapid transformation of all of the gas in a typi- 
cal galaxy disk into stars. Introducing a kinetic component 
to the feedback prolongs the depletion timescales some- 
what, but in all cases, best results are obtained by choos- 
ing long star formation timescales, i.e., values of c* ~ 0.05. 
For such low values of c*, e v must be less than 0.2 in order 
to prevent slowing down star formation to rates signifi- 
cantly lower than observed in Kennicutt 's empirical corre- 
lations. Similar constraints on e v can be derived from high 
resolution ID simulations of supernova remnants (Thorn- 
ton et al. 1998). We shall hereafter refer to the (rather 
extreme) choice of = 0.05 and e v — 0.2 as the "ki- 
netic feedback" case. Adopting c* > 0.05 and e v < 0.2 



produces results that are intermediate between the "mini- 
mal" and "kinetic" feedback cases we report here. Tests on 
isolated disks of varying circular velocity indicate that the 
"kinetic feedback" adopted here reduces substantially the 
efficiency of star formation in systems with circular veloc- 
ity below ~ 100 km s _1 , but has a more modest influence 
on the star formation rates in more massive systems, in 
rough agreement with current interpretation of observa- 
tional data (Martin 1999). 

3.3. The Initial Conditions 

We investigate two variants of the Cold Dark Matter 
scenario. The first is the former "standard" CDM model, 
with cosmological parameters ft = 1, h = 0.5, A = 0, 
normalized so that at z = the rms amplitude of mass 
fluctuations in 8ft -1 Mpc spheres is as — 0.63. Although 
this model fails to reproduce a number of key observations, 
such as the CMB fluctuations detected by COBE, it re- 
mains popular as a well-specified cosmological testbed and 
as a well studied example of a hierarchical clustering model 
of galaxy formation. The second is the currently popular 
low-density CDM model that includes a non-zero cosmo- 
logical constant and which is normalized to match COBE 
constraints and current estimates of the Hubble constant, 
£!o = 0.3, A = 0.7, h = 0.7, og = 1-1- Both models assume 
a value of ilb = 0.0125 h~ 2 for the baryon density pa- 
rameter of the universe, consistent with constraints from 
Big-Bang nucleosynthesis of the light elements (Schramm 
& Turner 1997). 

In both cosmologies, we simulate regions that evolve to 
form dark halos with circular velocities in the range (80, 
350) km s _1 at z — 0. These regions are selected from 
cosmological simulations of large periodic boxes and are 
rcsimulated individually including the full tidal field of 
the original calculation. All resimulations have typically 
32, 000 gas particles and the same amount of dark matter 
particles. The size of the resimulated region scales with 
the circular velocity of the selected halo, so that most sys- 
tems have similar numbers of particles at z = 0, regardless 
of circular velocity. This is important to ensure that nu- 
merical resolution is approximately uniform across all sys- 
tems. Gas particle masses range from 2.5 x 10 6 /i _1 M to 
9 x 10 7 h~ 1 M Q , depending on the system being simulated. 
Dark matter particle masses are a factor (fio — Qb)/&b more 
massive. Their mass is low enough to prevent artificial sup- 
pression of cooling due to collisional effects (Steinmetz & 
White 1997). Runs start at z = 21 and use gravitational 
softenings that range between 0.5 and 1.0 h^ 1 kpc. 

We have concentrated our numerical efforts on the 
sCDM scenario: about 60 galaxy models satisfy the con- 
ditions for analysis outlined below (§3.4). For compari- 
son, only seven ACDM galaxy models are considered here. 
This is because of our realization during the course of this 
study that the cosmological scenario has very modest ef- 
fects on the scaling laws we discuss here. Indeed, regarding 
disk scaling laws, galaxy models formed in the sCDM and 
ACDM scenarios differ mainly because of the adoption of 
different values of the Hubble constant. We discuss this in 
more detail below. 

3.4. Identification and Analysis of Model Galaxies 

Model galaxies are easily identified in our runs as star 
and gas "clumps" with very high density contrast. We 
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retain for analysis only galaxies in halos represented with 
more than 500 dark particles inside the virial radius; most 
of these systems have more than 1,000 star particles in 
each galaxy. This list is culled to remove obvious ongo- 
ing mergers and satellites orbiting the main central galaxy 
of each halo. The properties of the luminous component 
are computed within a fiducial radius, r ga i = 15 (Va/220 
km s _1 ) h^ 1 kpc. This radius contains all of the baryonic 
material associated with the galaxy and is much larger 
than the spatial resolution of the simulations. The rota- 
tion speeds we quote are also measured at that radius, 
although in practice the circular velocity in the models 
is rather insensitive to the radius where it is measured: 
similar results are obtained using r ga i = 3.5 (Va/220 km 
s" 1 ) h- 1 kpc (see Figure 1 of SN1). 

Although the resolution of the modeled galaxies is ade- 
quate to compute reliably quantities such as the total mass 
or circular velocity as a function of radius, it is still insuffi- 
cient to gain insight into the morphology of the galaxy. As 
a result, our sample is likely to contain a mixture of Hub- 
ble types, from SO to Sd (most of the simulated galaxies 
retained for analysis are largely rotationally supported). 

Galaxy luminosities are computed by simply adding the 
luminosities of each star particle, taking into account the 
time of creation of each particle and using the latest ver- 
sion of the spectrophotometric models of Bruzual & Char- 
lot (G.Bruzual & S. Chariot 1996, unpublished), see Con- 
tardo, Steinmetz & Fritze-von Alvensleben (1998) for de- 
tails. Corrections due to internal absorption and inclina- 
tion are neglected. The IMF is assumed in all cases to be 
independent of time. 
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Fig. 5. — Correlations between disk mass, Mdiak, disk rota- 
tion speed (measured at r ga i), Kot, absolute I-band magnitude 
at z = 0, Mi, and halo circular velocity, Va, found in the nu- 
merical experiments. Symbols are as in Figure 1. Solid line 
in the lower-left panel is the best fit to the observational data 
shown in Figure 1. Solid and dashed lines in the upper left 
panel correspond to constant disk mass-to-light ratio, Ti = 2.5. 
Solid line in upper-right panel outline the loci of galaxies that 
have assembled all available baryons. Slopes quoted in all pan- 
els correspond to unweighted least-squares fits. The rms scatter 
values correspond to the y-axis quantity relative to the least- 
squares fit, except for the Mi-V ro t panel, where it corresponds 
to the scatter in Mi. 



4. NUMERICAL SCALING LAWS 
4.1. The numerical Tully-Fisher relation 



The symbols with horizontal error bars in Figure 1 show 
the numerical Tully-Fisher relation obtained in our simu- 
lations. Solid squares and open circles denote the lumi- 
nosities and rotation speeds (measured at r ga {) of galaxy 
models formed in the sCDM and ACDM scenarios, respec- 
tively, under the "minimal feedback" assumption. Starred 
symbols correspond to the "kinetic feedback" case applied 
to the sCDM model. Error bars span the range in lumi- 
nosities corresponding to assuming either a Salpeter or a 
Scalo stellar initial mass function. 

A few points are clear from Figure 1. (i) The scat- 
ter and slope of the numerical Tully-Fisher relation are in 
good agreement with observation, (ii) The zero point of 
the numerical relation is offset by almost 2 (1.25) magni- 
tudes for sCDM (ACDM) galaxy models, (iii) The effects 
of kinetic feedback on the numerical Tully-Fisher relation 
are quite modest. Indeed, only a slight dimming in galax- 
ies with V TO t 100 km s^ 1 is clearly noticeable in Figure 
1. We analyze these three results in detail below. 

4.1.1. The Slope 

As discussed in §2.1, agreement between the slopes of 
the numerical and observed /-band Tully-Fisher relation 
imply a tight relation between the fraction of the mass as- 
sembled into the galaxy, / m dsk, the disk mass-to-light ra- 
tio, Tj, and the ratio between disk and halo circular veloc- 
ities, V ro t/VA (see eq. 5). Correlations between these pa- 
rameters obtained in our numerical simulations are shown 
in Figure 5. 

The bottom left panel of Figure 5 is identical to Fig- 
ure 1, but labels includes the (unweighted) best-fit slope 
and rms deviation (in magnitudes) of each numerical rela- 
tion. Symbols are as in Figure 1. As mentioned above, the 
slopes of the sCDM relations are consistent with the ob- 
served one. The same applies to the ACDM best-fit slope, 
which is virtually indistinguishable from the sCDM one. 
This result is only weakly affected by the narrow dynamic 
range covered by the ACDM simulations. Indeed, extend- 
ing our analysis to more poorly resolved clumps in the 
ACDM runs (systems with V Iot < 200 km s -1 , see dotted 
open circles in this panel) shows that there is little signifi- 
cant difference between the effective slopes obtained in the 
sCDM and ACDM models (see also Figure 2 of NS2). The 
good agreement between observed and numerical slopes in- 
dicate that /mdsk, Y_r, and Vrot/VA approximately follow 
the constraint prescribed by eq. 5. The rest of the panels 
in Figure 5 illustrate the nature of this relation. 

The upper left panel in Figure 5 shows that, to a large 
extent, the stellar mass-to-light ratios of galaxy models 
are approximately constant; the dotted and dashed lines 
in this panel correspond to T/ = 2.5 ^constant. On the 
other hand, the disk mass fraction is far from constant 
from system to system. As illustrated by the upper right 
panel in Figure 5, baryons in low circular velocity halos 
are much more effective at condensing into central galax- 
ies than those in more massive systems. Indeed, as Va 
increases, the numerical results move away from the solid 
line representing galaxies which have assembled all of the 
available baryons. 6 The fraction of baryons assembled 
into the central galaxy varies from ~ 70% in halos with 



6 The symbols corresponding to the ACDM model (open circles) have been rescaled downwards in the upper right panel of Figure 5 so that 
the Afjisk = ATJJY s °lid line is the same as in sCDM. Vertical deviations from this solid line thus indicate in both cases variations in the 
fraction of baryons assembled into the central galaxy. 
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y A <100 km s _1 to less than ~ 20% in V A ~ 300 km s" 1 
halos. This result is unlikely to be an artifact of limited 
numerical resolution since, as discussed in §3.3, the numer- 
ical resolution of the simulations do not depend system- 
atically on halo circular velocity. Furthermore, the same 
trend was found in the convergence study of Navarro & 
Steinmetz (NS97, see their Table 1); the observed trend 
between / m dsk and Va must therefore reflect the decreas- 
ing efficiency of gas cooling in halos of increasing mass 
discussed in §2.1.1. 

In spite of the large systematic variations in the disk 
mass fraction put in evidence by Figure 5, the numerical 
slope of the Tully-Fishcr relation is in good accord with 
observations because / m( jsk variations are compensated by 
corresponding changes in the V^/Va ratio. Low mass ha- 
los have higher / m dsk values, but also higher V^/Va ratios 
than more massive ones (bottom right panel of Figure 5). 
For constant T/, agreement with the observed TF slope 
requires / mdsk = M disk /M A oc M disk /Vl oc (V^/Va) 3 

(see eq. 5), so that if Aidisk oc V^, then V mt oc V^ 3 must 
be satisfied. Inspection of the exponents ( "slopes" ) listed 
in the right-hand panels of Figure 5 demonstrate that this 
condition is approximately satisfied in the numerical ex- 
periments. 



SC[)U: 6, -00 n.-OOO 
aCDM: e,=00. ^. = 005 
w:DM:f v =0.0, c.=0.0S 




Fig. 6. — As in Figure 2, but including the results of the 
numerical simulations. The dotted line shows the proportion- 
ality /mdsk oc (K-ot/VA) 3 - The thick dashed line shows the 
"adiabatic contraction" prediction assuming that a / m dsk oc 
V^ ' 7 , and that the NFW concentration parameter is given by 
c = 20(Va/100 km s _1 )~ 1//3 , as found in the numerical experi- 
ments. See text for further details. 



Figure 6 illustrates this conclusion more explicitly: 
the numerical results follow very closely the / m dsk oc 
(V^/Va) 3 proportionality outlined by the dotted line. 
As discussed in §2.1.1, this proportionality results from 
the combination of three effects: (i) the decreasing disk 
mass fraction in halos of increasing mass (/mdsk oc 
approximately, see upper-right panel of Figure 5), (ii) 
the decrease in concentration of halos of increasing mass 
(c « 20 (Fa/100 km s" 1 )- 1 / 3 , see NFW96 and NS2), and 
(hi) the response of the dark halo to the assembly of the 
disk. The thick dashed line in Figure 6 shows the / m dsk vs. 
(K-ot/VA) relation that results from these three premises 
using the adiabatic contraction approximation described 
in §2.1.1; it clearly reproduces the numerical results re- 
markably well. 

We conclude that the agreement between the slopes of 
the numerical and observed Tully-Fisher relations lends 
support to the view that the Tully-Fisher slope is the (non- 
trivial) result of the combined effects of the halo structure 



and of its response to the assembly of the baryonic com- 
ponent of the galaxy. 



4.1.2. The Zero Point 

In contrast with the good agreement found between 
models and observations for the Tully-Fisher slope, it is 
clear from Figure 1 that there is a serious mismatch in the 
zero-point of the numerical and observed Tully-Fishcr re- 
lations. This is reflected in Figure 6 as a systematic offset 
between the numerical data points and the thick curves 
labeled "sCDM" and "ACDM" (analogous to the lines in 
Figure 2, but drawn here for T/ = 2.5, as appropriate for 
our numerical results). At given / m dsk, the ratio V xo t/VA 
is much larger than required by observations, leading to 
galaxies which, at fixed rotation speed, are too faint to be 
consistent with observations, as clearly shown in Figure 1. 
Simulated galaxies formed in the sCDM and ACDM sce- 
narios are physically very similar, since they have similar 
disk mass-to-light ratios and share the same location in the 
/disk vs. (V^/Va) plane. The reason why ACDM mod- 
els appear to be in better agreement with observations in 
Figure 1 is largely due to the adoption of a higher value of 
Hubble's constant for this model. This issue is discussed 
in detail by NS2. From the analysis in that paper and the 
discussion in §2.1.1, it is clear that the zero-point discrep- 
ancy is testimony to the large central concentrations of 
dark halos formed in the two cosmologies we explore here. 



4.1.3. The Scatter 

According to eq. 4, the scatter in the numerical Tully- 
Fishcr relation is determined by the variance in the mass- 
to-light ratio, Tj, the disk mass fraction, / m dsk) and the 
velocity ratio, V ro t/V A . These can be derived from the 
rms values quoted in the right-hand panels and upper-left 
panel of Figure 5, respectively, and would imply, taken at 
face value, a large dispersion (>0.6 mag rms) for the nu- 
merical Tully-Fisher relation. This is actually about three 
times larger than measured in the simulations (see rms val- 
ues quoted in bottom left panel of Figure 5). The reason 
behind the small scatter is once again linked to the dynam- 
ical response of the halo to the disk assembly, as discussed 
in §2.1.1. At fixed halo mass (i.e., fixed Va), galaxies 
where fewer than average baryons collect into the central 
galaxy have lower than average V^/Va ratios, and vicev- 
ersa. This is shown in Figure 7, where we plot the residuals 
from the best-fit power laws to the data presented in the 
right-hand panels of Figure 5. As anticipated in §2.1.1, the 
residuals scale in such a way that variations in the frac- 
tion of baryons assembled into galaxies scatter along the 
observed Tully-Fisher relation, reducing substantially the 
resulting scatter. This helps to reconcile the large scatter 
predicted by analytical estimates (Eisenstein & Loeb 1996) 
with the results of our numerical simulations. Again, the 
halo structure and response to the assembly of the galaxy 
are crucial for explaining the observed properties of the 
Tully-Fisher relation. 
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Fig. 7. — Residuals from least squares fits to the numerical 
data presented in the right-hand panels of Figure 5. Solid line 
corresponds to the condition S log Mdisk = 35 log Va required 
so that galaxy models scatter along the Tully-Fisher relation, 
reducing substantially its scatter. 

4.2. The angular momentum of simulated disks 

In agreement with prior work (see, e.g., NS97 and ref- 
erences therein), we find that the baryonic components of 
the simulated galaxy models are quite deficient in angular 
momentum. This is easily appreciated in Figure 3, which 
shows that, at a given V TO t, the specific angular momen- 
tum of observed disks exceeds that of numerical models 
by more than one order of magnitude. This is due to the 
transfer of angular momentum from the baryons to the 
halo associated with merger events during the formation 
of the disk, as first suggested by Navarro & Benz (1991), 
and later confirmed by Navarro, Frenk & White (1995) and 
NS97. Indeed, as shown in Figure 4, the baryonic compo- 
nents of simulated galaxies have retained, on average, less 
than 15% of the specific angular momentum of their sur- 
rounding halos, placing them in the (fbdsk, fj) plane well 
outside of the constraints imposed by the observed scaling 
laws between luminosity, rotation speed, and disk size. 

4.3. The effects of feedback 

Our implementation of feedback seems to have only a 
very modest impact on the results discussed above, even 
for the rather extreme "kinetic feedback" case we tried. 
In the case of the Tully-Fisher relation only systems with 
Vrot ^ 100 km s _1 are affected (see, e.g., starred symbols in 
Figure 1), and the net effect is actually to make the zero- 
point disagreement worse: the rotation speeds of these low- 
mass galaxies are roughly unchanged, but feedback makes 
them fainter. This shows yet again that the zero-point 
problem afflicting the Tully-Fisher relation is not a result 
of processes that affect the baryonic component, such as 
feedback, but of the high central concentration of dark 
matter near the centers of dark halos and the relatively 
high disk mass-to-light ratios of our models. 

Feedback also seems to have a relatively minor impact 
on the "angular momentum problem" as well. Disks sim- 
ulated with our "minimum" and "kinetic" feedback im- 
plementations have angular momenta well below those of 
observed spirals (Figures 3 and 4). Feedback seems able to 
slow down the transformation of gas into stars in low-mass 
systems but fails to prevent most of the gas from collaps- 
ing into dense disks at early times. The final galaxy is 
thus built up as the outcome of a hierarchical sequence 
of merger events during which the gas transfers most of 
its angular momentum to its surrounding halo (NS97 and 
references therein). 



We conclude that the feedback algorithm explored here 
is unable to bring the angular momenta of simulated galax- 
ies into agreement with observations. Gas must be effec- 
tively prevented from collapsing into, or removed from, 
dense disks at early times, a requirement that requires a 
much more efficient transformation of supernova energy 
into gas bulk motions than our algorithm can accomplish 
(Weil, Eke & Efstathiou 1998). However, it is difficult to 
see how this can be attained without violating constraints 
posed by observed correlations between gas density and 
star formation rates in isolated disks (Kennicutt 1998). As 
discussed in §3.2, the "kinetic feedback" model parameters 
adopted here are already quite extreme; for example, in- 
creasing e v beyond our choice of 20% in an effort to expel 
more gas from star forming regions would lead to mod- 
els that disagree substantially with Kennicutt 's findings 
(SN3). Although it cannot be discounted that feedback 
implementations different from the one adopted here may 
lead to improved results, we conclude that accounting si- 
multaneously for the luminosity, velocity, and angular mo- 
mentum of spiral galaxies in hierarchical models remains 
an unsolved problem for the CDM cosmogonies we explore 
here. 

5. SUMMARY AND DISCUSSION 

We report here the results of numerical experiments de- 
signed to explore whether cosmologically induced correla- 
tions linking the structural parameters of cold dark mat- 
ter halos may serve to elucidate the origin of disk galaxy 
scaling laws. The numerical experiments include gravity, 
pressure gradients, hydrodynamical shocks, radiative cool- 
ing, heating by a UV background, and a simple recipe for 
star formation. Feedback from evolving stars is also taken 
into account by injecting energy into gas that surrounds 
regions of recent star formation. A large fraction of the 
energy input is in the form of heat, and is radiated away 
quickly by the dense, cool, star-forming gas. The remain- 
der is introduced through an extra acceleration term that 
affects directly the kinematics of the gas in the vicinity of 
star forming regions. 

We find that the slope of the numerical Tully-Fisher re- 
lation is in good agreement with observation, although not, 
as proposed in previous work, as a direct result of the cos- 
mological equivalence between mass and circular velocity 
of dark halos. Rather, the agreement results from a deli- 
cate balance between the dark halo structure, the fraction 
of baryons collected into each galaxy, and the dynamical 
response of the dark halo to the assembly of the luminous 
component. 

Massive halos are significantly less efficient at assem- 
bling baryons into galaxies than their low-mass counter- 
parts, a trend that agrees with expectations from theo- 
retical models where the mass of the central galaxy is de- 
termined by the efficiency of gas cooling. As a result of 
the structural similarity of dark halos, systems that col- 
lect a large fraction of available baryons into a central 
galaxy have their rotation speeds increased substantially 
over and above the circular velocity of their surrounding 
halos. The combined effect leads, for Cold Dark Matter 
halos, to a direct scaling between disk mass and rotation 
speed, Mdisk oc V^ ot , which, for approximately constant 
stellar /-band mass-to-light ratio, is in very good agree- 
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mcnt with the observed slope of the /-band Tully-Fishcr 
relation. 

The scatter in the numerical Tully-Fisher relation (~ 
0.25 mag rms, see lower left panel in Figure 5) is sub- 
stantially smaller than observed, a somewhat surprising 
result given the sizeable dispersions in disk mass-to-light 
ratios and disk mass fractions, as well as in the ratios be- 
tween disk rotation speed and halo circular velocity ob- 
served in our numerical simulations. The small scatter is 
mainly due to the fact that, although disk mass fractions 
may vary widely from galaxy to galaxy, these variations 
are strongly correlated with variations in the disk rotation 
speed. Model galaxies therefore scatter along the Tully- 
Fishcr relation, minimizing the resulting scatter (Figure 
7). This is again a non-trivial result stemming from the 
detailed dark halo structure and from its response to the 
assembly of the luminous component of the galaxy. 

Docs the agreement extend to passbands other than 
17 Since our galaxy models have approximately constant 
mass-to-light ratio in the /-band, the simulated slope will 
be approximately similar in all passbands where mass ap- 
proximately traces light, such as in the infrared K band. 
This is consistent with the results of Verheijen (1997), who 
find a slope of —7.0 in the /-band and of —7.9 in K, when 
his complete sample of galaxies is analyzed. The A'-band 
slope is significantly steeper than the /-band's (—10.4 in 
K versus —8.7 in /) only when a restricted sample with 
favorable inclination, smooth morphology, steep HI profile 
edges, and free from bars is considered. The strong depen- 
dence on sample selection reflects the large uncertainties 
in the numerical values of the slope that result from the 
small dynamic range in velocity (less than a factor of ~ 3 
typically) covered by existing Tully-Fishcr samples. Tak- 
ing this into account, and considering that the selection 
criteria of Verheijen's complete sample is closer to those 
in our analysis (§3.4), we conclude that our models are 
generally consistent with the slopes of the / and K band 
Tully-Fisher relations. 

The Tully-Fisher relation is also known to be signifi- 
cantly shallower in bluer passbands (see, e.g., Verheijen 
1997) a result that can be traced to systematic variations 
in the mass-to-light ratios as a function of disk rotation 
speed. Qualitatively our simulations can reproduce the 
trend, although it is quantitatively much weaker than ob- 
served (SN1). This is because fewer stars are formed late 
in our models than implied by observations (SN3). This 
effect may have led us to underestimate the scatter in Tj. 
However, considering that the measured scatter is much 
smaller than the observed one, we still consider our con- 
clusion that slope and scatter are consistent with observa- 
tions to be very solid. 

In contrast with this success, we find that the zero-point 
of the numerical Tully-Fisher relation is offset by about 
~ 2 (~ 1.25) magnitudes relative to observations for galax- 
ies formed in the sCDM (ACDM) scenario. The offset can 
be traced to the high-central concentration of dark matter 
halos formed in these cosmogonies and, to a lesser extent, 
to the relatively high stellar mass-to-light ratios found in 
our numerical models (Tj *=s 2.5). The "concentrations" of 
dark matter halos would have to be reduced by a factor of 
~ 3-5 or, alternatively, T/ would have to be reduced to val- 
ues as low as ~ 0.4-0.5 in order to restore agreement with 
the observed zero-point. As discussed by NS2, neither al- 



ternative is quite palatable, since they involve substantial 
modifications either to the underlying cosmological model, 
to our understanding of the spectrophotometric evolution 
of stars, or to our assumptions about the initial stellar 
mass function. Furthermore, mass-to-light ratios as low 
as that would exacerbate the problem posed by the small 
fraction of mass and the large fraction of angular momen- 
tum collected by the central disk in low-density universes 
(see §2.2 and Figure 4). 

Including the energetic feedback from evolving stars and 
supernovae does not improve substantially the agreement 
between the properties of model galaxies and observation. 
This is because, in order to fit the empirical correlation 
linking star formation rate to gas surface density, our feed- 
back model is only efficient in low-mass systems with cir- 
cular velocities below ~ 100 km s _1 . The progenitors of 
massive spirals typically exceed the threshold circular ve- 
locity at high redshift, leading to efficient gas cooling and 
to early onset of star formation. The effects of the kinetic 
feedback implementation on the Tully-Fishcr relation are 
therefore minor, and only noticeable in systems with low 
circular velocities. The zero-point disagreement actually 
worsens at the faint end as star formation is slowed down, 
affecting the absolute magnitudes of model galaxies more 
than their rotation speeds. The limited impact of our feed- 
back algorithm also implies that the luminous component 
of model galaxies is assembled through a sequence of merg- 
ers, accompanied by a substantial loss of its angular mo- 
mentum to the surrounding dark matter. The angular 
momentum of model galaxies is, as a consequence, about 
one order of magnitude less that that of observed spirals, 
in agreement with previous work (cf. NS97). 

6. CONCLUDING REMARKS 

The results we discuss here illustrate the difficulties 
faced by hierarchical models that envision the formation of 
disk galaxies as the final outcome of a sequence of merger 
events. Agreement with observations appears to require 
two major modifications to our modeling, (i) Dark ha- 
los that are much less centrally concentrated than those 
formed in the two Cold Dark Matter scenarios we explore 
here, (ii) Feedback effects dramatically stronger than as- 
sumed here, affecting substantially the cooling, accretion, 
and star formation properties of dark halos perhaps as 
massive as Va ~ 200-300 km s^ 1 . 

The extreme feedback mechanism apparently required 
to bring the mass and angular momentum of disk galaxies 
in accordance with hierarchical galaxy formation models 
should have major implications on further observational 
clues and traces of the galaxy formation process. Given 
the paucity of present-day examples of this process at work 
(Martin 1999), we are led to speculate that star formation 
(and therefore feedback) episodes were much more violent 
in the past than in the local universe, perhaps as a re- 
sult of the lower angular momenta and increased surface 
density of star forming regions at high redshift. Large 
scale winds driven by early starbursts, perhaps associated 
with the formation of stellar spheroids, may rid sites of 
galaxy formation of early-accreting, low-angular momen- 
tum baryons altogether, allowing higher-than-average an- 
gular momentum material to collapse later into the mod- 
estly star-forming, extended disks prevalent in the local 
universe. 
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It remains to be seen whether such scenario for the as- 
sembly and transformation of baryons into galaxies with- 
stands observational scrutiny, but so far the z ~ 3 evi- 
dence being gathered from "Ly-break" galaxies seems to 
point to a past where starbursts may have been the norm 
rather than the exception (Heckmann 1999). Resolving 
the puzzle created by disk galaxy scaling laws has thus 
the potential to unravel questions of fundamental impor- 
tance to our understanding of the assembly and evolution 
of galaxies in a cosmological context. 
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